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A NEW  OPTIMALITY  CRITERION  METHOD  FOR 
LARGE  SCALE  STRUCTURES* 


M.R.  Khan*,  K.D.  Willmert**,  and  W.A.  Thornton*** 
Clarkson  College  of  Technology 
Potsdam,  New  York 


Abstract 

An  optimality  criterion  method,  which  exploits 
the  concept  of  one  most  critical  constraint,  is 
reported.  The  method  eliminates  the  need  to  calcu- 
late a large  set  of  Lagrange  multipliers  for  the 
active  constraints,  and  also  eliminates  the  need 
for  a decision  as  to  whether  or  not  a particular 
constraint  should  be  considered  active.  The  method 
can  treat  multiple  load  conditions  and  stress  and 
displacement  constraints.  Application  of  the 
method  to  a number  of  truss  and  frame  structures 
demonstrates  the  efficiency  and  accuracy  of  the 
method. 


formula  for  redesign.  If  there  is  only  one  type  of 
constraint  (ie,  either  stress  or  displacement  or 
buckling) , the  redesign  process  requires  only  an 
analysis  of  the  structure  and  an  application  of  the 
appropriate  recursion  formula.  There  is  no 
requirement,  in  addition  to  an  analysis  of  the 
structure,  to  solve,  (a)  a set  of  linear  algebraic 
equations  for  a set  of  Lagrange  multipliers  (as  in 
Ref.  8) , or  (b)  to  solve  a linear  program  based  on 
a linearization  of  an  assumed  set  of  active  and 
potentially  active  constraints  (as  in  Ref.  1) , or 
(c)  to  solve  a nonlinear  programming  problem  (Ref. 

3)  in  the  active  and  potentially  active  constraints 
in  order  to  obtain  a new  design. 


I.  Introduction 

The  problem  of  structural  optimization  has 
become  of  great  interest  to  many  researchers  during 
the  past  few  years.  The  goal  of  this  recent  work 
has  been  primarily  to  obtain  a minimum  weight  struc- 
ture subject  to  various  constraints  in  minimal  com- 
putational time  and  with  minimal  computer  storage. 
The  efficiency  of  earlier  painfully  slow  mathema- 
tical programming  techniques  for  large  structural 
problems  has  been  improved  considerably  by  Schmit, 
Farshi,  and  Miura  ' ' . Venkayva,  Gellatly,  Berke, 
Knot,  Gorzynski  and  Thorn  ton4 ' ^ ' 7 have  developed 

physical  optimality  criterion  techniques  to  effi- 
ciently design  large  scale  structures.  Also,  Dobbs 
and  Nelson,  and  Rizzi®'^  have  recently  used  mathe- 
matical optimality  criterion  methods  based  on  the 
Kuhn-Tucker  conditions  to  obtain  minimum  weight 
designs  efficiently.  Khan,  Thornton  and  Willmerti0 
applied  efficient  physical  optimality  criterion 
techniques  to  simple  structures  and  complex  high 
speed  mechanisms. 

The  development  of  the  method  presented  here 
was  motivated  by  a desire  to  extend  to  problems 
with  multiple  constraints  of  different  types  (ie, 
stress  and  displacement  constraints)  the  simplicity 
inherent  in  physical  optimality  criterion  methods 
developed  for  single  constraints  of  each  type.  For 
instance,  the  stress  ratio  method  has  over  the  years 
demonstrated  a remarkable  ability  to  efficiently 
produce  minimum  weight  designs  or  near  minimum 
weight  designs  for  a groat  variety  of  multiloaded 
structures  under  stress  constraints.  Likewise, 
physical  optimality  criterion  methods  for  displace- 
ment constraints,  have  been  derived  and  applied 
with  success. 

Each  of  these  independent  physical  optimality 
criterion  methods  gives  rise  to  a simple  recursion 
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In  this  paper,  recursion  formulas  for  stress 
and  displacement  constraints,  which  result  from 
the  Kuhn-Tucker  necessary  conditions  for  each  type 
of  constraint,  are  incorporated  into  a design 
algorithm  which  exploits  the  concept  of  a single 
most  critical  displacement  constraint.  The 
algorithm  requires  only  one  analysis  of  the  struc- 
ture per  design  cycle.  Redesign  of  each  member  is 
achieved  by  means  of  one  of  two  recursion  formulas. 
No  sets  of  Lagrange  multipliers  need  be  calculated, 
no  subsidiary  LP  or  NLP  must  be  solved,  no  decision 
as  to  active  or  potentially  active  constraints  must 
be  made,  and  no  move  limits  need  be  used.  The 
method  is  applicable  to  two  and  three  dimensional 
trusses  and  two  dimensional  frames,  of  fixed 
geometry,  under  multiple  load  conditions  and  stress 
and  displacement  constraints. 


II.  Theory 

The  design  problem  to  be  solved  here  can  be 
stated  as:  find  the  vector  of  design  variables 
A = (A^ , A2»  ...»  A^|)  such  that  the  volume  of  the 
structure 
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where  A^  and  1^  are  the  cross-sectional  area  and 
length  of  the  ith  member,  N is  the  number  of  mem- 
bers, 0^  is  the  stress  in  the  ith  member  in  the 
kth  load  condition,  K is  the  number  of  load  condi- 
tions, and  0|  is  the  limiting  stress  in  the  ith 
member.  Also,  ujfc  is  the  displacement  in  jth  con- 
strained degree  of  freedom,  u j is  the  limiting 
value  of  the  displacement  in  the  jth  constrained 
degree  of  freedom,  and  J is  the  number  of  displace- 
ment constrained  degrees  of  freedom. 
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Stress  Constraints 


Combining  eqs.  (8)  and  (9) , results  in 


Considering  stress  constraints  alone,  the 
Kuhn-Tucker  conditions  for  the  design  problem  of 
eqs.  (1)  and  (2)  results  in  the  well  known  stress 
ratio  formula  for  redesign  (see  Ref.  10  for 
example) 
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where  v is  the  iteration  counter.  If  design 
variable  linking  is  used  to  form  groups  of  design 
variables,  where  members  of  one  group  are  the  same 
size,  eq.  (3)  is  applied  to  each  member  of  a group, 
and  the  largest  from  eq.  (3)  is  taken  as  the 
size  for  all  members  of  the  group. 


Displacement  Constraints 

Considering  displacement  constraints  alone, 
the  Lagrangian  for  the  design  problem  of  eqs.  (1) 
and  (2)  is 

J K 

L * V ♦ l l A (u  - u > (4) 

j = l k-1  3*  3 3 

and  the  Kuhn-Tucker  necessary  conditions  for  a 
minimum  are 
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which  is  the  optimality  criterion  which  must  be 
satisfied  at  the  optimum  design.  From  eq.  (10) 
the  following  recursive  formula  results 
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If  design  variable  linking  is  used  and  n members 
are  to  have  the  same  design  variable  A^,  eq.  (11) 
is  written  as 


**  a.  y i. 
1 i=i  3 


(12) 


where  the  summation  over  j in  eq.  (12)  is  over 
those  members  which  have  common  design  variable  A^ . 

In  eqs.  (11)  and  (12),  n is  a relaxation 
parameter  which  is  used  to  control  the  rate  of 
convergence  and  stability  of  the  method.  It  is  the 
only  arbitrary  parameter  involved  in  the  algorithm. 
Values  between  0.001  and  0.2  have  been  found  to  be 
appropriate . 
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Suppose  now  that  the  pth  constrained  displacement 
in  the  qth  load  condition  is  exactly  active,  and 
the  other  constrained  displacements  are  not.  Then 
eq.  (5)  becomes 
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By  means  of  the  unit  load  theorem  of  structural 
analysis,  the  derivative  in  the  first  of  eqs.  (6) 
can  be  written  (Ref.  11) 
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where  is  the  stiffness  matrix  of  the  ith  mem- 
ber, Xi  is  the  displacement  vector  for  the_^ith 
member  due  to  the  qth  load  condition,  and  x^  is 
the  displacement  vector  for  the  ith  member  due  to 
a unit  load  applied  at  the  location  and  in  the 
direction  of  the  pth  constrained  degree  of  free- 
dom. Substituting  eqs.  (1)  and  (7)  into  (6) 
gives: 

A.l.  - A xTK.x.  - 0 i*l N (8) 

i i pq  l i i 

and  sunning  eq.  (8)  over  all  members  results  in 


The  derivation  presented  above  is  incorporated 
into  the  design  algorithm  of  the  next  section.  Its 
use  is  justified  for  multi-displacement  constrained 
problems  because  of  the  selection  of  the  most 
active  (or  most  violated)  constraint.  In  this 
method,  as  well  as  most  currently  available  tech- 
niques, there  is  normally  only  one  most  active 
constraint  at  any  iteration.  There  may  be  many 
constraints  which  are  nearly  active — this  of  course 
is  especially  true  at  the  optimal  design,  but  only 
one  which  is  most  active.  The  true  optimal  design 
may  be  one  having  several  active  constraints,  but 
this  is  almost  never  exactly  obtained.  This 
characteristic  is  further  enhanced  by  the  fact  that 
finite  arithmetic  is  used,  so  round-off  eliminates 
additional  equalities.  In  the  special  case  where 
two  or  more  displacement  constraints  are  exactly 
equal  because  of  symmetry  or  other  structural  limi- 
tations, these  exactly  equal  displacements  are 
treated  as  one  constraint.  In  the  method  presented 
here,  this  most  active  (or  most  violated)  con- 
straint in  some  load  condition  is  considered  to  be 
the  only  active  constraint:  all  other  displacement 
constraints  are  considered  inactive. 

The  recursion  relations  of  eqs.  (11)  and  (12) 
have  been  applied  to  several  displacement  con- 
strained problems,  but  practical  problems  will  be 
those  with  both  stress  and  displacement  con- 
straints. Thus,  the  stress  recursion  formula  of 
eq.  (1)  has  been  combined  with  the  displacement 
recursion  formulas  of  eqs.  (11)  and  (12)  and  an 
important  scaling  procedure  to  produce  a design 
procedure  which  is  applicable  to  stress  and  dis- 
placement constrained  trusses  and  frames  under 
multiple  load  conditions. 


i 
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Also,  the  Lagrange  multiplier  associated  with 
the  critical  displacement  is  computed  from 
eq.  (9)  as: 


r 


III.  Design  Algorithm 

1.  Choose  any  uniform  design  Aa  , i-1,2 N. 

Choose  a value  of  the  relaxation  parameter  (r|)  . 
(Say  0.1S  to  0.08] 

2.  Analyze  the  design  for  each  load  condition. 

3.  Check  displacements  in  each  load  condition  at 
those  nodes  where  displacement  limitations  are 
imposed  and  determine  the  node  and  direction 
for  which  the  calculated  displacement  most 
closely  approaches  (or  exceeds)  the  allowable 
displacement.  This  is  the  most  critical  dis- 
placement (upq) • 

4.  Knowing  the  magnitude  of  the  most  critical  dis- 
placement (Upq)  from  step  3 and  the  value  of 
the  allowable  displacement  (Up) , scale  the 
chosen  design  so  that  the  most  critical  dis- 
placement becomes  active.  All  other  displace- 
ment constraints  will  then  be  inactive.  Let 
the  scaled  design  be  denoted  by  , where 

I u I 

A)  - — 23-  A i-1,2, ...,N  (13) 

1 |-  . l 


If  the  structure  was  analyzed  with  the  scaled 
design  then  displacement  vectors  calculated  at 
step  2 would  have  been 


i-1,2. .. . .N 


(14) 


and  stiffne?-  matrix  from  the  scaled  design 
would  be: 


(15) 


5. 


From  the  scaled  displacement  vectors  (Xj)  and 
design  A^,  compute  the  maximum  stress,  max|o.  | 

k ** 

in  each  member  i.  Also,  determine  the  stress 
response  ratio  for  each  member  and  let  the  most 
critical  response  ratio  be  obtained  for  the  nth 
member.  This  is  denoted  by  R„.  If  R„  > 1 
N 

compute  V,  - R ( ) All.),  or  if  R <1  compute 
1 n . i i n — 

„ l-l 


i-1 
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6.  Using  the  scaled  design,  apply  a unit  load  only 
at  the  node  and  in  the  direction  of  the  active 
displacement  constraint.  Let  the  set  of  ^ 
resulting  nodal  displacements  be  denoted  by  x^. 
Note  that  this  is  the  only  unit  load  that 
needs  to  be  applied,  and  that  the  structural 
stiffness  matrix  inverted  at  step  2 is  used 
here  as  scaled  in  step  4 to  compute  x^. 


From  eq.  (7)  compute 
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8.  Group  the  members  as  follows: 

1 • If  3Up^/3A|  < 0 or  CL  > 0^ , member  i 
belongs  to  group  Gj 

ii.  Otherwise,  member  i belongs  to  Group  G^ 

Note  that  either  group  could  be  empty  and 
a particular  member  would  belong  to  only  one 
group  at  a time. 


Use  the  stress  ratio  formula,  eq.  (3)  to 
resize  the  elements  of  G^,  as: 
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10.  Resize  the  elements  of  Gj  using  eq.  (11)  (or 
eq.  12)  , as 
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11.  Scale  the  design  and  compute  the  new  critical 
response  ratio  R^  and  new  using  ste^s  2 
through  5.  If  the  quantity  |(Vj-V^)/VjJ  is 
less  than  e ( a small  number  ranging  between 
0.001  to  0.010),  then  go  to  step  12:  other- 
wise check  the  following 

(a)  If  Vj  < Vi  continue  with  step  6 with  the 
old  value  of  0. 

(b)  If  Vj  > V1#  Rj,  < 1,  and  R^  > 1 , the 
designer  may  stop  at  this  point  and  the 
design  of  the  previous  iteration  would  be 
very  close  to  the  optimum  design.  Other- 
wise p is  reduced  to  one  third  or  one 
quarter  of  the  starting  value  and  the 
process  is  continued  with  step  6. 

12.  If  the  converged  design  of  step  5 is  com- 
pletely displacement  dominated  then  R^,  would 
be  less  than  1 and  this  design  is  the  optimum 
design.  If  the  converged  design  of  step  5 is 
completely  stress  dominated,  that  is,  all 
members  are  in  G^  and  hence  overstressed, 
simply  scale  the  design  (multiplying  all  the 
design  variables  with  R^)  so  that  no  stress 
constraint  is  violated  to  achieve  the  optimum 
design.  If  the  converged  design  has  some  mem- 
ber overstressed  while  others  not,  then 
following  situations  may  occut: 

(a)  If  Vj  < vi  and  (R„-l)  < 0.05,  scale  the 
design  of  the  previous  iteration  by  multi- 
plying all  the  design  variables  with  Rn 
and  this  is  then  taken  to  be  optimum 
design. 
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(b)  If  < VL  and  (R^-i)  < 0.05.  scAie  the 
design  of  the  current  iteration  by  multi- 
plying All  the  design  variables  by  R^  and 
this  is  then  taken  to  be  the  optimum 
design. 

(o  if  <r;-i>  > 0.05  reduce  the  value  of  n to 
half  or  one  third  of  the  starting  value 
and  go  back  to  step  6 and  repeat  the 
process. 

Choice  of  Relaxation  Parameter  (Q) 

This  is  the  only  arbitrary  parameter  in  the 
design  procedure.  It  controls  the  stability  and 
convergence  of  this  method.  Experience  indicates 
that  a value  of  n between  0.001  to  0.2  results  in 
optimum  designs  beinq  obtained  without  difficulty. 
It  is  important  to  note  that  selecting  the  value 
fran  this  range  does  not  affect  the  optimum  «iesign. 
The  same  design  will  be  obtained  using  any  value  of 
H between  0.001  and  0.2,  but  it  will  be  located  in 
fewer  iterations  with  the  larger  values.  One  dif- 
ficulty with  the  larger  values  of  f)  is  that  the 
technique  brings  the  design  close  to  the  optimum  in 
a very  few  analyses,  but  oscillations  will  occur 
very  close  to  the  optimum.  This  is  easily  detected 
when,  at  a particular  iteration,  the  scaled  design 
weighs  more  than  the  previous  design.  When  this 
occurs,  n is  reduced  and  the  procedure  is 
stabi li zed. 


IV.  Results 

In  this  section,  results  for  six  classical 
truss  examples  and  two  frame  examples  are  presented. 
These  are  intended  to  show  the  efficiency  and 
accuracy  of  the  design  algorithm  of  Section  III. 

1.  Ten  Bar  Truss 

This  is  a cantilever  truss  which  has  been 
studied  by  many  researchers  (Ref.  1, 2, 4, 5,8,9) . It 
is  shown  in  Pig.  1.  The  material  is  aluminum  of 
specific  weight  0 * 0.1  lb/in ^ and  modulus  of  elas- 
ticity E • 10xl06  psi.  Displacement  limits  of 
♦ 2.0  inches  are  imposed  on  all  nodes  in  both  direc- 
tions, and  the  limiting  value  of  stress  in  each  mem- 
ber is  ♦ 25,000  psi.  No  desiijn  variable  linking  is 
used  so  there  are  10  independent  design  variables. 
Two  cases  are  considered.  Case  1 has  * 100  K, 

* 0,  and  case  2 has  * 150  K,  P 2 3 50  K.  A 
single  loading  condition  is  considered  in  each  case. 
A lower  limit  on  member  size  of  0.1  in2  is  enforced. 

The  final  design  for  case  1 is  given  in  Table 
la.  In  this  case  the  problem  was  started  with  a 
uniform  design  with  each  cross-sectional  area  equal 
to  100  in2.  A starting  value  of  n 8 0.2  was  chosen 
and  was  automatically  changed  to  0.05  as  the 
design  came  close  to  the  optimum.  At  iteration  15 
a weight  of  5085  lbs  was  obtained  and  the  design 
was  similar  to  one  previously  reported  by  other 
researchers.  However,  the  algorithm  did  not  stop 
automatically  until  iteration  18  at  which  point  the 
weight  dropped  to  5067  lbs,  the  displacement  of 
node  1 in  the  y-direction  was  -2.0  inch,  the  dis- 
placement of  node  4 in  the  y-direction  was  0.4% 
below  its  limiting  value,  and  member  5 had  stress 
2.71%  below  its  yield  value.  It  is  interesting  to 
note  that  the  final  design  has  the  lowest  weight 
ever  achieved  for  this  problem. 


The  final  design  for  case  2 is  given  in  Table 
lb.  This  problem  was  started  with  the  same  initial 
design  and  n value  as  for  case  1.  The  design  was 
automatically  converged  at  iteration  9 when 
member  5 had  its  stress  equal  to  the  limiting  value 
and  the  displacement  of  node  4 in  the  y-direction 
was  0.3%  below  its  specified  limit.  The  final 
design  obtained  is  in  good  agreement  with  previous 
designs. 

2.  Four  Bar  Space  Truss 

The  structure  is  a four  bar  pyramid  truss 
shown  in  Fig.  2.  The  material  is  aluminum  with 
p 3 0.1  lb/in^  and  E * 10xl06  psi.  Stress  limits 
of  ♦ 25,000  psi  are  imposed  on  all  members.  No 
design  variable  linking  is  used.  TVo  cases  are 
considered.  Case  1 has  a loading  of  Px  3 10  K, 

Py  * 20  K,  and  Pz  * -60  K,  and  a displacement 
limit  of  t 0.3  inch  is  imposed  at  the  top  joint  in 
the  z-direction.  Case  2 has  a loading  of  Px  * 40  K, 
Py  * 100  K,  Pz  3 -30  K,  and  displacement  limits  at 
the  top  joint  are  ♦ 0.3  inch  in  the  x-direction, 
i 0.5  inch  in  the  y-direction  and  t 0.4  inch  in 
the  z-direction.  Results  are  qiven  in  Table  2. 

This  table  shows  good  correspondence,  with  previous 
results,  of  the  design  obtained  with  the  new  method, 
and  its  efficiency.  The  init'al  desiqn  for  both 
cases  had  all  members  at  100  in2.  In  both  cases  1 
and  2 member  3 had  stress  equal  to  its  limitinq 
value,  while  in  case  1,  displacement  of  the  top 
node  in  the  z-direction  was  3.fl%  below  its  limit 
and  in  case  2 displacement  in  the  y-direction  was 
1.9%  below. 

3_. Twenty-two  Member  Space  Truss 

This  structure,  which  is  shown  in  Fig.  3,  has 
each  joint  connected  to  every  other  joint  by  a 
member,  except  that  members  between  support  joints 
are  excluded.  It  was  studied  in  Ref.  12  in  the 
context  of  determining  the  global  optimum  of 
trusses  with  vanishing  members. 

All  members  are  aluminum  with  E = 10x10^  psi 
and  p * 0.1  lb/in^.  The  22  members  are  linked 
into  7 groups  as  shown  in  Table  3.  Table  3 also 
qives  the  limitinq  stresses  for  each  group  of  mem- 
bers. Displacement  constraints  of  t 0.2  inches  at 
all  nodes  in  all  directions  are  imposed,  and  a 
minimum  member  size  of  0.1  in2  holds  if  a member  is 
not  prescribed  to  vanish.  Three  load  conditions, 
as  given  in  Table  4,  are  considered  in  each  of  3 
design  cases.  Case  1 has  all  qroups  of  members 
nonvanishinq , case  2 has  the  members  of  group  4 set 
to  zero,  and  in  case  3,  the  members  of  group  3 
vanish.  Table  5 summarizes  the  results  of  the 
3 cases  obtained  by  the  method  of  this  paper  and 
compares  them  with  the  results  of  Ref.  12.  Case  1 
is  the  global  optimum  for  thii.  truss.  The  present 
method  achieves  a desiqn  with  weight  within  1%  of 
the  global  minimum  weight  in  5 analyses.  Cases  2 
and  3 converge  to  designs  very  close  to  the 
results  of  Ref.  12  in  just  6 analyses. 

The  initial  design  for  all  3 cases  was  uniform 
with  all  members  at  10  in2.  The  starting  values  of 
parameter  n for  the  three  different  cases  were 
arbitrarily  chosen  to  be  0.2,  0.125  and  0.1  respec- 
tively. These  changed  to  one  quarter  of  their 
startinq  values  at  the  end  of  optimization 
process.  Also, the  desiqn  process  was  studied 
by  starting  all  3 cases  with  the  same  value  of  n. 


The  final  designs  obtained  were  the  same  as  those 
presented  in  Table  3. 

4 . Twenty- five  Bar  Transmission  Tower  Truss 

This  much  studied  truss  (Refs.  1.2,4, 5, 8, 9)  is 
shown  in  Fig.  4.  The  material  of  all  members  is 
again  aluminum  with  E * 10x10”  psi  and  p • 

0.1  lb/inJ.  Design  variable  linking  is  used  to 
reduce  the  number  of  independent  design  variables 
from  25  to  8.  Table  6 gives  the  members  of  each 
design  variable  group.  CHiis  problem  was  solved 
using  both  25  and  8 independent  design  variables, 
with  insignificant  differences  in  CPU  time.  The 
results  for  the  8 design  variable  case  are  pre- 
sented here  for  purposes  of  comparison  with  pre- 
vious results.)  The  stress  limits  for  each  group 
of  members  are  also  given  in  Table  6.  Displacement 
limits  of  4 0.35  inch  are  imposed  on  every  node  in 
every  direction.  TVo  load  conditions  are  con- 
sidered. These  are  given  in  Table  7.  Table  8 
gives  the  final  design  obtained  and  compares  this 
with  pteviously  obtained  designs.  The  comparison 
indicates  that  the  method  of  this  paper  gave  a 
design  similar  to  those  previously  obtained,  but 
with  a weight  about  2%  higher.  The  problem  was 
started  with  H equal  to  0.1  and  all  members  at  100 
in^.  The  design  automatically  converged  at  9 itera- 
tions with  horizontal  displacements  at  the  joints  1 
and  2 equal  to  their  limiting  values.  The  final 
design  is  completely  displacement  dominated. 

5 . Seventy-two  Hember  Space  Truss 

This  structure,  shown  in  Fig.  5,  has  been 
studied  previously  in  Kefs.  1,2, 4, 5 and  6.  All 
members  are  aluminum  with  F.  * 10x10^  psi  and  p - 
0.1  lb/in3.  Stress  limits  of  4 25,000  psi  are 
imposed  on  all  members.  Displacement  limits  of 
4 0.25  inch  in  the  x and  y directions  are  imposed 
on  the  4 top  nodes.  A lower  limit  of  0.1  in^  is 
imjosed  on  all  members.  Design  variable  linking  is 
used.  Members  are  placed  in  16  qroups  as  shown  in 
T.ible  9.  Thus,  there  are  16  independent  design 
variables.  Two  load  conditions  are  considered. 

These  are  given  in  Table  10.  Table  11  gives  the 
final  designs  obtained  for  two  initial  values  of  n, 
and  compares  these  with  previous  results.  The 
design  procedure  was  started  with  all  members  equal 
to  100  in^.  Starting  with  n " 0.15,  it  was  noted 
that  at  iteration  8 a weight  of  394  lbs  was 
achieved  but  the  procedure  continued  until  itera- 
tion 10  when  it  was  automatically  stopped  with  a 
weight  of  388  lbs.  At  the  optimum,  in  the  second 
load  condition  the  first  four  members  had  their 
stress  equal  to  their  limiting  values  while  the 
displacements  of  node  1 in  the  x and  y directions 
were  2.1%  below  their  specified  limits. 

6.  TVo  Hundred  Member  Planar  Truss 

This  structure,  previously  studied  in  Ref.  13, 
is  shown  in  Fig.  6.  All  members  are  steel  with  E * 
JOxlO6  ps i and  P * 0.283  lb/in3.  A stress  limit  of 
4 10,000  psi  is  imposed  on  all  members,  and  dis- 
placement limits  of  4 0.5  inch  are  imposed  on  all 
nodes  in  both  directions.  The  structure  is  symme- 
tric about  the  vertical  centerline.  This  reduces 
the  number  of  independent  design  variables  to  105. 
Three  load  conditions  are  considered: 

1.  1 K in  positive  y direction  at  all  nodes 
on  line  AB; 


2.  10  K in  negative  z direction  at  all  nodes 
on  lines  AB,  CD,  EF,  GH,  and  IJ; 

3.  load  conditions  1 and  2 acting  together. 

The  final  design  obtained  is  given  in  Table  12. 

The  final  weight  of  32,996  lbs  obtained  with  8 
analyses  and  34.35  minutes  of  CPU  time  on  an  IBM 
360/65  compares  favorably  with  the  weight  of 
31,020  lbs  obtained  in  90  minutes  of  CPU  time  on  an 
IBM-7094- II-7044-DCS . Comparing  the  design 
obtained  by  the  present  method  with  that  obtained 
by  Ref.  13  indicates  that  they  are  somewhat  dif- 
ferent. Results  of  several  solutions  obtained  by 
the  method  of  this  paper  indicate  that  the  region 
of  the  optimum  is  flat,  i.e.,  designs  of  signifi- 
cantly varying  member  sizes  are  possible  for 
essentially  the  same  weight. 

Both  designs  of  Table  12  have  one  displacement 
constraint  active  at  the  optimum.  This  is  the 
displacement  at  node  I in  the  z-direction. 

7.  Three  Member  Frame 

The  structure  is  shown  in  Fig.  7.  It  is  a 
three  member  rigid  frame.  Each  member  is  treated 
by  one  finite  element.  Axial,  shear,  and  bending 
moment,  are  included  in  the  formulation,  resulting 
in  6 degrees  of  freedom  per  element  and  3 degrees 
of  freedom  per  joint.  The  material  is  steel  with 
E * 30x10^  psi  and  p = 0.283  lb/in^.  The  design 
variable  for  each  member  is  the  cross-sectional 
area  A.  The  section  modulus  S and  moment  of. 
inertia  I are  related  to  area  as  S = 9A  and  I = 

75A.  These  relationships  were  chosen  to  give 
sections  representative  of  available  wide  flange 
shapes  while  maintaining  the  linearity  among  A,  S, 
and  1.  The  stress  limits  for  all  members  are 
4 24,000  psi.  One  load  condition,  as  shown  in  Fig. 

7,  is  imposed.  Three  cases  are  considered.  Cases 
1 and  2 include  the  above  stress  limits  and  the 
following  displacement  constraints;  case  1 has  the 
displacements  of  joints  2 and  3 limited  to  4 0.2 
inch  in  the  x and  y directions  and  case  2 has  the 
same  di splacements  limited  to  ± 0.07  inch.-  For 
case  3,  the  stress  limits  are  ignored  and  only  dis- 
placement constraints  of  4 0.2  inches  at  joints  2 
and  3 in  both  directions  are  considered.  Table  13 
gives  the  results  of  these  3 cases  and  compares 
them  to  previously  obtained  results.  It  can  be 
seen  that  excellent  agreement  has  been  obtained  at 
a fraction  of  the  CPU  time  required  for  these  pre- 
vious results. 

Initial  designs  for  Briggs  (Ref.  14)  and  SUMT 
were  uniform  at  75  in^,‘  and  those  for  the  method  of 
this  paper  uniform  at  100  in^.  The  n values  of 
Table  13  were  constant  during  the  design  process. 

8.  Twenty-five  Member  Frame 

The  structure  is  shown  in  Fig.  8.  Members  are 
defined  as  in  Example  7.  One  load  condition  is 
considered  as  shown  in  Fiq.  8,  and  one  finite  ele- 
ment is  used  per  member.  All  members  are  100 
inches  in  length  except  the  diaqonal  members  which 
are  141.4  inches  long.  Stress  limits  are  4 24,000 
psi  for  all  members.  Two  cases  are  considered. 

Case  1 has  the  above  stress  limits  and  displace- 
ment limits  of  4 3.0  inches  at  joints  1,  2,  3,  4, 

5,  and  6 in  both  directions.  Case  2 has  the  above 
stress  limits  and  displacement  limits  of 
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♦ 0.05  inch  at  joints  1,  2,  3,  4,  5 and  6 in  both 
directions.  The  minimum  member  size  is  5 in2. 
Results  are  shown  in  Table  14.  Both  cases  were 
started  with  all  members  equal  to  100  in2.  Case  1 
is  compared  with  results  from  Ref.  14  with  excel- 
lent agreement  in  the  designs.  The  method  of  this 
paper  can  be  seen  to  produce  the  optimal  design 
with  a drastic  reduction  in  the  CPU  time  required 
for  the  method  of  Ref.  14.  The  design  for  case  1 
is  fully  stressed  at  the  optimum  and  the  displace- 
ment limits  are  inactive.  The  case  2 design  is 
displacement  constrained,  with  no  active  stress 
constraints.  No  previous  results  were  available 
for  comparison. 

The  n values  given  in  Table  14  did  not  change 
during  the  design  process. 

IV.  Conclusions 

A new  design  algorithm  has  been  developed  for 
stress  and  displacement  constrained  trusses  and 
frames  under  multiple  loadings.  By  means  of  an 
extensive  set  of  test  problems,  the  method  as  been 
shown  to  be  both  accurate  and  efficient.  In  all 
problems  studied,  known  results  were  reproduced 
very  closely  with  the  number  of  structural  analyses 
required  in  the  iterative  process  approximately  the 
same  as  the  number  required  by  the  current  most 
efficient  methods.  When  it  is  considered  that  the 
computational  effort  required  per  iteration  for  the 
method  of  this  paper  is  considerably  less  than  that 
required  for  all  other  current  methods,  and  also 
that  the  core  storage  required  is  essentially  only 
that  required  for  the  analysis  capability,  the 
present  method  can  be  seen  to  be  very  simple  as 
well  as  being  highly  efficient. 
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Figure  3.  Twenty-two  Member  Space  Truss 


Table  la.  Comparison  of  Final  Designs  for  Ten  Bar  Truss,  Case  1 


Final  Cross-Sectional 

Areas  (in2) 

Member 

Schmit 

i Miura 

Schmit  & 

Venkayya 

Gellatly 

Dobbs  t 

Rizzi 

This 

No. 

NEWSUMT 

CONMIN 

Farshi 

fc  Berke 

Ne 1 son 

Paper 

Ref.  2 

Ref.  2 

Ref.  1 

Ref. 4 

Ref. 5 

Ref . 0 

Ref. 9 

i 

30.670 

30.57 

33.432 

30.416 

31.350 

30.500 

30.731 

30.980 

2 

0.100 

0.  369 

0. 100 

0.128 

0.100 

0.100 

0. 10 

0.10 

3 

23.760 

23.97 

24.260 

23.408 

20.030 

23.290 

23.934 

24.163 

4 

14.590 

14.73 

14.260 

14.904 

15.600 

15.428 

14.733 

14.805 

5 

0.100 

0. 10 

0.100 

0.101 

0.140 

0.100 

0.100 

0.100 

6 

0.100 

0.  364 

0.100 

0.101 

0.240 

0.210 

0.100 

0.406 

7 

8.570 

8.547 

8.  338 

8.696 

8.350 

7.649 

8.542 

7.547 

8 

21.070 

21.11 

20.740 

21 .084 

22.210 

20.980 

20.954 

21.046 

9 

20.960 

20.77 

19.690 

21.077 

22.060 

21.818 

20.036 

20.937 

10 

0.100 

0.320 

0.100 

0.186 

0.100 

0.100 

0.100 

0.100 

Wt(lbs) 

5076.85 

5107.3 

5089.0 

5004.9 

5112.0 

5080.0 

5076.66 

5066.98 

Analyses 

13 

14 

24 

26 

19 

15 

11 

184 

aA  weight  of  5085  lbs  was  achieved  after  15  analyses 


Table  lb.  Comparison  of  Final  Designs  for  Ten  Bar  Truss,  Case  2 


Final  Cross-Sectional 

Areas  (in^) 

Member 

Schmit 

& Miura 

Schmit  & 

Venkayya 

Oel latly 

Dobbs  & 

Rizzi 

This 

No. 

NEWSUMT 

Ref  .2 

CONMIN 

Ref  .2 

Farshi 

Ref.l 

Ref. 4 

& Berke 

Ref  .5 

Ne 1 son 

Ref  .8 

Ref  .9 

Paper 

i 

23.550 

23.55 

24.289 

25.190 

. 

25.813 

23.533 

24.716 

2 

0.100 

0.176 

0.100 

0.  363 

- 

0.100 

0.100 

0.100 

3 

25.290 

25.20 

23.346 

25.419 

- 

27.233 

25.291 

26.541 

4 

14.360 

14.39 

13.654 

14.327 

- 

16.653 

14.374 

13.219 

5 

0.100 

0.100 

0. 100 

0.417 

- 

0.100 

0.100 

0. 108 

6 

1.970 

1.967 

1.969 

3.144 

- 

2.024 

1.9697 

4.835 

7 

12. 390 

12.400 

12.670 

12.083 

- 

12.776 

12.389 

12.664 

8 

12.810 

12.860 

12.544 

14.612 

- 

14.218 

12.825 

13.775 

9 

20.340 

20.410 

21.971 

20.261 

22.137 

20.328 

18.438 

10 

0.100 

0. 100 

0.100 

0.513 

0.100 

0.100 

0.10 

Wt(lbs) 

4676.96 

4684.11 

4691.84 

4895 . 60 

- 

5059.7 

4676.92 

4792.52 

Analyses 

11 

10 

23 

13 

* 

12 

12 

9 

Table  2a 

Final  Designs,  Four  Bar  Pyramid,  Case  1 


Member 

No. 

Final  Cross 

-Sectional  Areas  (in2) 

Schmit  6 
Farshi 

Ref . 1 

Venkayya 

Ref  .4 

Th  is 

Paper 

i 

0.0 

0.277 

0.0 

2 

3.765 

4.1527 

3.651 

3 

0.769 

0.746 

0.769 

4 

2.514 

2.477 

2.759 

Wt(lbs) 

117.89 

126.43 

121.50 

Analyses 

16 

37 

6 

Table  2b 

Final  Designs,  Four  Bar  Pyramid,  Case  2 


Member 

No. 

Final  Cross-Sectional  Areas  (in2) 

Schmit  & 
Farshi 

Ref  .1 

Venkayya 

Ref. 4 

This 

Paper 

i 

3.210 

3.147 

3.419 

2 

2.614 

2.147 

2.511 

3 

2.159 

2.162 

2.159 

4 

0.0 

0.0 

0.0 

Wt(lbs) 

128.53 

128.561 

130.625 

Analyses 

14 

- 

7 

Table  3 

Member  Linking  Groups  and  Stress  Limits, 
Twenty- two  Member  Space  Truss 


Table  4 

Load  Conditions  for  Twenty- two 
Member  Space  Truss 


Design 

Variable 

Group 

Number 

Members  of 
Group 

Lower 

Limiting 

Stress 

(psi) 

Upper 

Limiting 

Stress 

(psi) 

1 

1,2. 3, 4 

24,000 

36, 

000 

2 

5,6 

30,000 

3 

7,8 

28,000 

4 

9,10 

26,000 

s 

11,12,13,14 

22,000 

/ 

6 

15,16,17,18 

20,000 

7 



19,20,21,22 

18,000 

36, 

000 

Lend 

Load  Components 

Condition  Node 

Px 

py 

Number 

(Kips) 

(Kips) 

(Kips) 

i i 

El 

0 

-5 

2 

Bi 

0 

-5 

3 

-20 

0 

-30 

mM' 

-20 

0 

-30 

2 1 

-20 

-5 

2 

-20 

-50 

3 

-20 

-5 

4 

-20 

-50 

3 1 

-20 

0 

35 

2 

-20 

0 

0 

3 

-20 

0 

0 

4 

-20 

0 

-35 

Table  5.  Final  Design  Comparison,  Twenty-two  Member  Space  Truss 


Group 

Case 

i 

Case 

2 

Case  3 

Number 

Sheu  & Schmit 

This  Paper 

Sheu  t Schmit 

This  Paper 

Sheu  & Schmit 

This  Paper 

Ref.  12 

Ref.  12 

Ref.  12 

i 

2.6288 

2.5627 

2.6101 

2.5262 

2.5657 

2.4902 

2 

1.1624 

1.5530 

1.4234 

1.9529 

1.1331 

1.8126 

3 

0. 3433 

0.2813 

0.587 

0.5475 

0.0 

0.0 

4 

0.4231 

0.5124 

0.0 

0.0 

0.6461 

0.6581 

5 

2.7823 

2.6261 

2.7861 

2 . 5900 

2.6738 

2.5442 

6 

2.1726 

2.1314 

2.0891 

2.2178 

2.1768 

2.2419 

7 

1.9523 

2.2128 

2.0935 

2.2630 

2.1613 

2.2799 

Wt(lbs) 

1024.80 

1034.74 

1028.07 

1040.51 

1029.35 

1040.47 

Analyses 

.a 

5 

_a 

6 

_a  ^ 

6 

Not  Applj  cable 


Table  6 

Member  Linking  Groups  and  Stress  Limits, 
Twenty-five  Member  Transmission  Tower  Truss 


Table  7 

Load  Conditions,  Twenty-five 
Member  Transmission  Tower  Truss 


Design 

Variable 

Group 

Number 

Members  of 
Group 

Lower 

Limiting 

Stress 

(lbs/in2) 

Upper 

Limiting 

Stress 

(lbs/in2) 

1 

i 

35092.0 

40,000.0 

2 

2, 3, 4, 5 

11590.0 

3 

6, 7, 8, 9 

17305.0 

4 

10,11 

35092.0 

5 

12,13 

35092.0 

6 

14,15,16,17 

6759.0 

7 

18,19,20,21 

6959.0 

N 

f 

8 

22,23,24,25 

11082.0 

40,000.0 

Load 

Node 

Direction 

Condition 

X 

y 

z 

1 

i 

1 K 

10  K 

-5  K 

2 

0 

10  K 

-5  K 

3 

.5  K 

0 

0 

6 

.5  K 

0 

0 

2 

1 

0 

20  K 

-5  K 

2 

0 

-20  K 

-5  K 

• -■  - . 


55 


Table  8.  Final  Designs,  T\/enty- fi ve  Member  Transmission  Tower  Truss 


Members 
In  Group 
No. 


Analyses 

Needed 


Schmit 

& Miura 

NEWSUMT 
Ref  .2 

CON  MIN 
Ref  .2 

0.010 

1.985 

2.996 

0.010 

0.010 

0.684 

1.677 

2.662 

0.  166 
2.017 

3.026 

0.087 

0.097 

0.675 

1.636 

2.669 

545.172 

548.475 

10 

9 

Schrait  & 
Farshi 
Ref.l 


Final  Cross  Sectional  Areas  (in-*) 


Venkayv-  | Gellatly 


Dobbs  & 
Nelson 
Ref  .8 


Gellatly 
& Berke 
Ref.S 


0.0100 

2.0069 

2.9631 

0.0100 

0.0100 

0.6876 

1.6784 

2.6638 


545.36  553.4 


0.01 
1.9884 
2.9914 
0.01 
0.01 
0.684 
1 .6767 
2.6627 


545.163 


Areas  not  reported 

Table  9 

Member  Linking  Groups, 

Seventy- two  Member  Truss 

Table  10 
load  Conditions, 
Seventy- two  Member  Truss 

Design  Variable 
Group  Number 


Members  in  Group 


Load 

Condition  I Node 


Direction 


1,2, 3,4 

5,6,7,8,9,10,11,12 

13,14,15,16 

17,18 

19,20,21,22 

23,24,25,26,27,28,29,30 

31,  32,33,  34 

35,36 

37, 38, 39,40 

41, 42, 43, 44, 4S, 46, 47, 48 

49,50,51,52 

53,54 

55,56,57,58 

59,60,61,62,63,64,65,66 

67,68,69,70 

71,72 


Table  11.  Final  Designs,  Seventy-two  Member  Truss 


Members  of 


Final  Cross-Sectional  Areas  (in^) 


Group 

Schmi t t 
NEWSUMT 
Ref  .2 

. Miura 

CONMIN 
Ref. 2 

Schmi t L 
Farshi 

Ref  .1 

Venkayya 

Ref  .4 

Gellatly 
l Berke 
Ref  .5 

Berke 

6 Knot 
Ref  .6 

This 

Paper 

(n-o.i) 

This 

Paper 

(0-0.15) 

1 

0.1565 

0.1558 

0.1585 

0.161 

0.1492 

0.1571 

0.1494 

0.1519 

2 

0.5458 

0.5484 

0.5936 

0.557 

0.7733 

0.5385 

0.5698 

0.5614 

3 

0.4105 

0.4105 

0. 3414 

0.  377 

0.4534 

0.4156 

0.4434 

0.4378 

4 

0.5699 

0.5614 

0.6076 

0.506 

0. 3417 

0.5510 

0.5192 

0.5317 

5 

0.5233 

0.5228 

0.2643 

0.611 

0.5521 

0.5082 

0.6234 

0.5814 

6 

0.5173 

0.5161 

0.5480 

0.532 

0.6084 

0.5196 

0.5231 

0.5273 

7 

0.1000 

0.1000 

0.1000 

0.100 

0.1000 

0.1000 

0.100 

0.100 

8 

0.1000 

0.1133 

0.1509 

0.100 

0.1000 

0.1000 

0.1963 

0.1583 

9 

1.267 

1.268 

1.1067 

1.246 

1.0235 

1.2793 

1.2076 

1.2526 

10 

0.5118 

0.5111 

0.5792 

0.524 

0.5421 

0.5149 

0.5208 

0.5244 

11 

0.1000 

0.1000 

0.1000 

0.100 

0.1000 

0.1000 

0.100 

0.100 

12 

0.1000 

0.1000 

0.1000 

0.100 

0.1000 

0.1000 

0.100 

0.100 

13 

1.885 

1.885 

2.0784 

1.818 

1.4636 

1.8931 

1.7927 

1.8589 

14 

0.5125 

0.5118 

0.5034 

0.524 

0.5207 

0.5171 

0.5223 

0.5259 

15 

0.1000 

0.1000 

0.1000 

0.100 

0.1000 

0.1000 

0.100 

0.100 

16 

0.1000 

0.1000 

0.1000 

0.100 

0.1000 

0.1000 

0.100 

0.100 

Final 

1 

Wt(lbs) 

379.640 

379. 792 

388.63 

381.2 

395.97 

379.67 

386.718 

387.67 

Analyses 

Needed 

9 

8 

22 

12 

9 

5 

13 

10 

Final  Cross- 
Sectional 
Area  (in^) 


Venkayya 


Member 

Number 


Briggs 
Ref. 14 


19.74 

105.38 

30.13 


Final  Cross 
Sectional 
Area  (in^) 


Table  14 

Final  Design  Comparison  for 
Twenty-five  Member  Frame 


Venkayya 
Ref.  13 


4.798 
4.798 
5.662 
5.662 
5.737 
5.737 
6.688 
6.688 
6.274 
6.274 
7.285 
7.285 
5.695 
5.695 
6.713 
6.713 
8.989 
8.989 
20.687 
20.687 
9.594 
9.594 
1.156 
2.278 
3.346 
4.495 
5.626 
6.  770 
7.822 
8.969 
I 9.800 
10.95 


8 

574 

8 

574 

7 

864 

7 

864 

8 

40 

8 

40 

6 

545 

6 

545 

7 

062 

7 

062 

8 

095 

8 

095 

20 

046 

20 

046 

9 

454 

9 

454 

1 

860 

2 

397 

3 

762 

4 

191 

5 

799 

6 

252 

7 

107 

7 

520 

8 

038 

7 

913 

Member  Cross-Sectional  Areas 

tin2) 

Case  1 

Case  2 

Member 

Briggs 

This  Paper 

This  Paper 

Number 

Ref. 14 

(0=0.1) 

(0=0.1) 

1 

138.00 

129.55 

337.79 

2 

148. 58 

153.26 

293.07 

3 

154.08 

151.29 

162.15 

4 

28.34 

31.63 

69.68 

5 

128.93 

133.64 

170.02 

6 

5.00 

5.00 

52.19 

7 

130.10 

131.58 

217.43 

8 

15.72 

23.37 

108.73 

9 

162.97 

170.80 

233.83 

10 

5.00 

5.00 

110.55 

11 

120.97 

119.91 

170.29 

12 

111.06 

110.78 

181.92 

13 

5.00 

5.00 

37.35 

14 

122.00 

123.06 

109.72 

15 

5.00 

5.00 

105.77 

16 

52.96 

54.00 

147.31 

17 

191.76 

190.80 

233.56 

18 

5.00 

5.00 

191.63 

19 

119.70 

123.13 

336.97 

20 

5.00 

5.00 

199.56 

21 

123.74 

119.  32 

465.20 

22 

8.61 

5.00 

191.92 

23 

5.00 

5.00 

88.26 

24 

5.00 

5.00 

84.14 

25 

48.78 

48.67 

95.59 

boh 

187421 

188215 

463523 

Analyses 

_a 

15 

10 

IK3&M 

1849.00 

69.02 

38.79 

Not  Applicable 

b 

All  tines  on  IBM  360/65 


Table  13.  Final  Design  Comparison  for  Three  Member  Frame 


19.68 

105.43 

30.12 


This  Paper 
(0=0.15) 


19.81 

105.39 

30.18 


SUMT  This  Paper 
(0=0.15) 


17.78 

130.07 

69.31 


This  Paper 
(0=0.2) 


Analyses 


Not  Applicable 
bAll  times  on  IBM  360/65 


